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Abstract. This paper deals with the convergence/divergence issue of the Chapman-Enskog series expansion of the shear 

Oh 1 and normal stresses for a granular gas of inelastic hard spheres. From the exact solution of a simple kinetic model in the 

■^T uniform shear and longitudinal flows, it is shown that (except in the elastic limit) both series converge and their respective 

radii of convergence increase with inelasticity. This paradoxical result can be understood in terms of the time evolution of the 
>0 , Knudsen number and the existence of a nonequilibrium steady state. 
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Introduction. In a Newtonian fluid the irreversible parts of the stress tensor Py are linear in the velocity gradients 
Vjiij. For instance, if the flow is incompressible (i.e., V • u = 0), P„ = — r]odu x /dy, where 770 is the shear viscosity 
and it has been assumed that duy/dx = 0. Analogously, in a compressible flow characterized by a rate-of-strain tensor 
Vjiij = (du x /dx)SixSj X one has P xx = p — (4/3)Tlodu x /dx, where p is the hydrostatic pressure and the bulk viscosity 
has been neglected (as happens in the low-density limit). In the case of a dilute gas, Newton's law is derived from 
the Boltzmann equation under the assumption that the characteristic hydrodynamic length (L) associated with the 
gradients of density (n), temperature (T), or flow velocity (u) is much larger than the mean free path I of the gas 
particles, i.e., £ <C L. More in general, the Chapman-Enskog (CE) method expresses the solution of the Boltzmann 
equation as an expansion in powers of the Knudsen number |ju| = £/L UJ]. The leading terms in the CE expansion 
yield the Navier-Stokes constitutive equations and provide explicit expressions for the transport coefficients (like the 
shear viscosity rjo). 

An important basic question is the nature (convergent versus divergent) of the CE expansion. To isolate the problem, 
\Q \ let us consider the following subclasses of the full CE series: 
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P xy = -£r lk (du x /dy) 2k+ \ P xx = p-\^' k {du x /dx) k+ \ (1) 

k=0 J k=0 
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where 770 = T]q is the shear viscosity, r][ is a Burnett coefficient, T71 and r\' 2 are super-Burnett coefficients, and so 
on. The two full CE series of P xa (with a — y and a — x, respectively) reduce to the partial series (Q} if (and only 
if) V,Mj = Y xa 8i a 8j x and Vn = VT = Vy xa — 0, i.e., the only nonzero hydrodynamic gradient is a uniform shear or 
longitudinal rate fxa = du x /da. These are precisely the conditions characterizing two well-defined physical states: 
the uniform shear flow (USF, a — y) [2|, |3|, |4j and the uniform longitudinal flow (ULF, a = x) [5, 6]. In both states 
the characteristic hydrodynamic length is L ~ y/2T /m/\y xa \ (where m is the mass of a particle and the temperature is 
measured in energy units), so that the Knudsen number becomes the absolute value of ji = jm/v, where v ~ \JTT /m/t 
is a characteristic collision frequency. Thus, the CE expansions ([TJ can be recast into the dimensionless forms 
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Px,/p = -tlF xy (n), P xx /p = l--^F xx { i i), F .(m)=E c ^ 2 *, Mm)=L<V, (2) 
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where c* = ('Hk/p)v 2k+l an d c 'k = ( 1 7fc/p) v • Despite the simple definitions of the USF and the ULF, no exact 
solution of the nonlinear Boltzmann equation is known for those states. However, the problem becomes solvable in the 
framework of the Bhatnagar-Gross-Krook (BGK) model kinetic equation yi |6(] and the solutions show that, except 
for Maxwell molecules, the CE expansions (Q]) are divergent Jj,|a]. 

So far we have implicitly assumed conventional gases made of particles that collide elastically. What about the case 
of granular gases? A granular gas is a large collection of macroscopic particles which collide inelastically and are 
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maintained in a fluidized state. It can be conveniently modeled as a gas of smooth inelastic hard spheres characterized 
by a constant coefficient of normal restitution a < 1 [7]. For such a model, the CE method has been applied to the 
inelastic Boltzmann equation and the Navier-Stokes transport coefficients have been derived [8]. Taking into account 
that the nonlinear subclasses (Q]i of the full CE expansions, as said before, do not converge when the gas is made of 
elastic hard spheres (a = 1) Jj,|a], and considering the inherently non-Newtonian nature of the steady USF of granular 
gases [9], it seems intuitive to expect that the CE series (fTJ are also divergent in the case of inelastic hard spheres. 
However, it turns out that, paradoxically, the series ((TJ do converge in the inelastic case, the radii of convergence 
increasing with increasing inelasticity tilftlllll . 

Rheological properties. In the USF, the mass and momentum balance equations yield constant density n and shear 
rate j xy . However, n(t) = n(0)/[l + jxx(0)t] and fxx(t) = 7xv(0)/[l + j xx (0)t] in the ULF, so that j xx > corresponds 
to an expansion of the gas, while jxx < corresponds to a condensation phenomenon yj, la]. In either case the ratio 
jxa{t)/n(t) is constant. As for the energy balance equation of inelastic hard spheres in the USF and ULF, one has 

d t T(t) = -(2%,/3n)P„(t) - C(t)T(t), (3) 

where C, is the so-called cooling rate 0], which is approximately given by £ = A(l — oc 2 )vq, Vq « n\fT being an 
effective collision frequency for elastic spheres. The cooling term on the right-hand side of Eq. ([3]) competes with the 
viscous heating term in the USF (jxyPxy < 0) and in the ULF with jxx < (note that P xx > 0), so that, depending on the 
initial state, the temperature either grows or decreases with time until a steady state is eventually reached [9, 11, 12], 
except in the elastic case (£ = 0). On the other hand, in the ULF with j^ > both terms on the right-hand side of Eq. 
(01 represent cooling effects, so T(t ) — > and no steady state exists. 

In order to analyze in detail the CE series (Q]), it is convenient to consider the following BGK-like kinetic model 
H13I1 of the inelastic Boltzmann equation: 

(d,+v-V)/=-v(/-/o) + (£/2)d v -[(v-u)/], (4) 

where / is the velocity distribution function, /o is the local version of the homogeneous cooling state distribution [TJ], 
Y^-Vq is an effective collision frequency, so that e = £/v = l( 



and v = ^y^Vq is an effective collision frequency, so that e = £/v = 1(1 — a). Taking moments in Eq. (0|i one gets 



USF 4S = vJiv;0^' ULF:^ = v,- ( v + C + 37^, (5) 

where p = nT. Equations ([3]) and (O constitute a closed set of equations for the evolution of {7\.Pr, ,P yy } (USF) or 
{T,P X x} (ULF). To describe the non-Newtonian hydrodynamic regime we must focus on the nonlinear dependence 
of the scaled quantities ^(^(f)) = ~Pxy{t) / 'p{t)jl(t) and F xx (n(t)) = ~(3/4)[P xx (t)/p(t) - l]/ju(f) as functions of 
jl{t) = j xa (t)/v(t) °< 1/ y/TJt). Elimination of time in favor of jJ. in Eqs. (0 and © yields a single second-order 
ordinary differential equation for F xy (n) [ 10] and a single first-order ordinary differential equation for F xx (jx) llllll . 

Results and discussion. The numerical solutions for a gas of elastic hard spheres (a = 1) and for a gas of highly 
inelastic hard spheres (a — 0.5) are plotted in Fig. QI a )- The arrows on the curves indicate the direction followed 
by the time evolution of jJ.(t). In the inelastic case (a < 1, £ ^ 0) the evolution leads to the steady-state point 
H s = ±v/3e72(l + e),F xy (n s ) = l/(l+e) 2 for the USF and & = -3e(l + e)/2(l+3e),F xx (^) = (1 +3e)/(l +e) 2 
for the ULF (provided that j xx < 0; otherwise, jj.(t) — > °°). On the other hand, in the elastic case (a = 1, S = 0) the 
temperature monotonically increases with time in the USF and in the ULF with j xx < and thus the Knudsen number 
asymptotically vanishes, i.e., ji(t) — > 0. The dependence of ji s on a is shown in Figure[Hb). 

According to the preceding discussion, the zero Knudsen number value {ji — 0) is a "repeller" of the time evolution 
of ji(t) both for USF and ULF in the inelastic case. For elastic collisions, however, the state ji — is an "attractor" 
of ji(t) (except in the ULF with j xx > 0). Expressed in other terms, the fixed point ji = is unstable against any USF 
or ULF perturbation, no matter how weak it is, if the particles are inelastic. On the other hand, ji = is stable 
for elastic collisions in the USF and in the ULF with j xx < 0. Since the CE expansions ([TJ and (O are carried 
out about (and measure the departure from) the reference homogeneous state (jx = 0), it follows that the CE series 
diverge if \\m t ^ m ji(t) = (elastic case), whereas they converge if fi(t ) flees from ji — (inelastic case). In the latter 
situation, in addition, the radius of convergence must coincide with the steady-state value ji s . Consequently, one must 
have lim^oc |q/c£ + i | = jti s 2 and lim^,*, \c' k /c' k+l \ ~ |/x v |. This expectation is confirmed by an exact evaluation of the 
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FIGURE 1. (a) Plot of F xy (fl) (solid lines) and F xx (fi) (dotted lines) for a — 1 and a = 0.5; the open circles denote the steady- 
state points, (b) Dependence of the steady-state value fl s on the coefficient of restitution a. (c) Ratios [ct/cfc+il and |c^/cj. + 1 | for 

a = 0.5; the horizontal dotted lines represent /lj and \fi s \, respectively. 

coefficients cj and c' k for k < 450, as illustrated by Fig. [He) for a = 0.5. In mathematical terms, the nonlinear viscosity 
functions F m (}i) possess a singularity located at ji = in the elastic case and at fi = fi s ^ in the inelastic case. 

To conclude, the paradoxical regularization by inelasticity of the CE series (Q]) is directly related to the time evolution 
of the Knudsen number |/i(f)|, and hence of the temperature, Eq. (0, and thus it is not an artifact of the kinetic model 
©. Also, the convergence or divergence of the series (Q]) does not depend on whether the system actually is in the 
USF, the ULF, or in any other state. The advantage of the USF and ULF is that the full CE series of the shear and 
normal stresses reduce to the partial series ([T), thus allowing one to explore their character in a detailed way. 
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